Negative binomial multiplicity distribution in proton-proton collisions in limited 
pseudorapidity intervals at LHC up to y/s= 7 TeV and the clan model 
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Experiments at the Large Hadron Collider (LHC) have measured multiplicity distributions in 
proton-proton collisions at a new domain of center-of-mass energy (y^) in limited pseudorapidity 
intervals. We analyze multiplicity distribution data of proton-proton collisions at LHC energies as 
measured by the Compact Muon Solenoid (CMS) experiment in terms of characteristic parameters 
of the Negative Binomial Distribution (NBD) function that has played a significant role in describ- 
ing multiplicity distribution data of particle production in high energy physics experiments, in the 
pre-LHC energy-range, in various kinds of collisions for a wide range of collision energy and for 
different kinematic ranges. Beside a single NBD, we apply the formalism of weighted superposition 
of two NBDs to examine if the multiplicity distribution data of CMS could be better explained. 
The weighted superposition of two NBDs indeed explain the distribution data better at the highest 
available LHC energy and in large interval of phase space. The two-NBD formalism further reveals 
that the energy invariance of the multiplicity distribution of the "soft" component of particle pro- 
duction in hadronic collisions is valid at LHC also, as it is at RHIC and Tevatron. We analyze the 
data further in terms of clan parameters in the framework of the two-NBD model. 

PACS numbers: 13.85.Hd 



I. INTRODUCTION 

Experimental study of multiparticle production in 
high-energy hadronic (proton-proton, pp or proton- 
antiproton, pp) collisions has reached a new high, in 
terms of energy, at the Large Hadron Collider (LHC) 
A fast growth in energy of collisions could be pos- 
sible due to significant advancement of collider tech- 
nology in the last few decades. From tens of GeV 
at Intersecting Storage Ring (ISR), hundreds of 
GeV [I at Super Proton Synchrotron (SPS), both at 
CERN, finally the center-of-mass energy (y/s) of colli- 
sions has reached thousands of GeV first at Tevatron 
at Fermilab Q and then at LHC at CERN [5-9]. Re- 
markably, for this wide range of collision energy, the 
two-parameter Negative Binomial Distribution (NBD) 
function, as given below in Eq. - (1), played major 
role in describing multiplicity distributions of produced 
charged particles. 
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where (n) is the average multiplicity and the param- 
eter k is related to dispersion D, (Z?^ = (n^) — (n)^) 

by 
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The charged particle multiplicity distributions in pp 
collisions at the ISR energies and in pp collisions at 
^/i = 540 GeV at the SPS, fit with NBD function 
satisfactorily in the full pseudorapidity, "q (where rj — 
—ln[tan{9/2)] and 9 is the polar angle of the particle 
with respect to the counterclockwise beam direction) 
space as well as in limited pseudorapidity intervals. 
But, at Vi= 900 GeV SPS energy, a single NBD func- 
tion could describe the data only for small pseudorapid- 
ity intervals at the mid-rapidity region, while for larger 
intervals, where shoulder-like structure appeared in the 
multiplicity distribution, a single NBD function turned 
out to be inadequate. Appearance of sub-structures 
in multiplicity distributions at higher energies and in 
larger pseudorapidity intervals has been attributed [lol - 
[l3| to weighted superposition or convolution of more 
than one functions representing more than one source 
or process of particle productions. Such sub-structure 
in SPS data at y/s = 900 GeV and in Tevatron data 
at a/s =1.8 TeV could be well explained by weighted 
superposition of two NBD functions [ll| . The NBD is 
quite pertinent for pp collisions at energies available at 
the Large Hadron Collider (LHC) also, as has been re- 
ported p, @| first by 'A Large Ion Collider Experiment' 
(ALICE). 

The study of high-energy particle collisions in pre- 
LHC energy-range, where soft processes of particle pro- 
ductions dominate, barring application of perturbative 
quantum chromodynamics (pQCD), depends mostly 
on phenomenological models. Many of these models 
[l0l - l22l |. which deal with multiplicity distributions of 
produced particles, interpret matching of various data 
of multiplicity distributions involving NBD function 
well within respective framework. Sometimes, the k- 
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parameter in NBD has very different meanings in some 
of these so far successful approaches indicating that the 
very wide occurrences of NBD in high energy experi- 
ments is not yet a well understood phenomenon. In 
the given scenario, a detailed study of multiplicity dis- 
tribution data of pp collisions at new LHC energies in 
terms of NBD would be worth carrying out for bet- 
ter understanding of the role of NBD in multiparticle 
production. 



II. MULTIPLICITY DISTRIBUTIONS IN 
PROTON-PROTON COLLISIONS AT LHC 

At LHC, multiplicity distributions in proton-proton 
collisions at center-of-mass energies ^/s = 0.9, 2.36 and 
7 TeV have been measured by different experiments, in 
different kinematic ranges and for different classes of 
events. All these LHC-experiments find that the mean 
multiplicities at the new LHC energies {\/s= 2.36 and 7 
TcV) had been underestimated by the event generators 
(like PYTHIA, PHOJET etc.) in use. 

The ALICE has measured primary charged particles 
at -^/s — 0.9 and 2.36 TeV in the mid-T; region in three 
limited overlapping ?7-intervals \ri\ < rjc = 0.5, 1.0 and 
1.3 0, in non-single diffractive (NSD) inelastic proton- 
proton collisions. At a/s = 7 TeV, instead of NSD in- 
elastic events, ALICE analyzed 6] an event class requir- 
ing at least one charged particle in \r]\ < 1 and mea- 
sured multiplicity distribution in that 77-interval only. 
The distributions measured by ALICE in the three t]- 
intervals, \r]\ < ijc — 0.5, 1.0 and 1.3, at the two en- 
ergies, = 0.9 TeV, ^/s = 2.36 TeV have been re- 
ported [3 to match fairly well with NBD. The NBD 
fit to the distribution at \/s= 7 TeV, measured by AL- 
ICE, has been reported p to be slightly underestimat- 
ing the data at low multiplicity {n < 5) and slightly 
overestimating the data at high multiplicity (n > 55). 
The Compact Muon Solenoid (CMS) experiment has 
measured primary charged hadrons for all the three 
LHC energies in non-single diffractive (NSD) inelas- 
tic proton-proton collisions in the mid-77cm region in 
five overlapping 77-intervals jryj < rjc = 0.5, 1.0, 1.5, 
2.0 and 2.4 [3| around the center-of-mass pseudora- 
pidity (r/cm = 0). The CMS experiment did not fit 
the distributions with NBD or other distribution func- 
tions but reported a change of slope in P„ for n > 20 
in its largest ?7-interval, rjc < 2.4. This feature be- 
comes more pronounced with increasing ^/s. A Toroidal 
LHC Apparatus (ATLAS) experiment has measured Q 
charged particle multiplicities for different event classes 
characterized by different lower cuts on the number of 
charged particles {rich < 1; 2 and 6) in different kine- 
matic ranges (pr > 100 MeV, 500 MeV in \t]\ <2.5). 
In contrast to ALICE, CMS and ATLAS experiments, 
which measured multiplicity distributions in the cen- 
tral region in pseudorapidity, the Large Hadron Col- 



lider Beauty (LHCb) experiment at LHC has its detec- 
tor with geometrical acceptance limited in the forward 
region. The LHCb experiment has analyzed multiplic- 
ity distributions for hard interaction events (at least one 
long track with transverse momentum, pT > 1 GeV/c) 
from pp interactions at y/s= 7 TeV in non-overlapping 
pseudorapidity bins of width |7/| < 0.5 in the pseudora- 
pidity range 2.5 < 77 < 4.5. 

It is important to note, at this point, that a few of the 
phenomenological models have already been contrasted 
with the LHC data. To explain the appearance of sub- 
structure in the distribution, the set of CMS data of 
charged hadron multiplicity has been analyzed j21|] in 
the framework of Independent Pair Parton Interaction 
(IPPI) [l^l which shows that the number of soft pair 
parton interactions from colliding particles and so the 
density of the partonic medium is large for the LHC 
data and increases with energy. Within the framework 
of the IPPI model, the findings favor enlarged role of 
collective effects in pp collisions at LHC. Similar con- 
clusion is obtained [2l| from analysis of the data in 
terms of Quark Gluon String Model (QGSM) fH ^ 
that fits better to the data than the IPPI model. An- 
other attempt [2^ to describe the multiplicity distribu- 
tion data at the new LHC energies has been in terms 
of a two-component model in quantum statistical ap- 
proach [lo| which described hadronic collision data up 
to ^/s = 540 GeV. The model considers convolution of 
a NBD (with k—1) function and a Poisson Distribution 
(PD)function representing two components of source of 
particle productions: a thermally equilibrated chaotic 
source and another coherent source, respectively. The 
study revealed poor agreement between the data and 
the model. However, normalized moments of multiplic- 
ity distributions for the LHC data have been repro- 
duced ^] by a model considering measured distribu- 
tion as superposition of a PD describing particle emis- 
sion from one source and a NBD describing distribution 
of other sources. 



III. OBJECTIVE 

Our objective is a detailed study of multiplicity dis- 
tributions of proton-proton collisions at available LHC 
energies in limited pseudorapidity intervals in terms of 
behavior of the characteristic parameters of NBD with 
respect to changing width of pseudorapidity interval, 
|?7| and center-of-mass energy, ^/s of collision. 

Of the LHC data, we chose to analyze the data 
recorded by the CMS experiment 0, for the reason 
that, (a) contrary to the other LHC experiments, CMS 
experiment has measured multiplicity distributions for 
all the three LHC energies {^/s = 0.9, 2.36 and 7 
TeV), available so far, in similar pseudorapidity inter- 
vals rjc = 0.5 to rjc = 2.4 (five intervals) and for the 
same class (non-single-diffractive or NSD) of events, fa- 
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cilitating systematic study of dependence on 77^ and ^/s 
for same class of events with more data points, (b) 
detailed study of CMS data in terms of NBD func- 
tion is still absent, (c) the phenomenological studies 
[21I [25I [26j with published data of multiplicity distri- 
butions froia proton^ proton collisions at LHC energies 
have primarily dealt with the CMS data and analysis of 
same data by different phenomenological models help 
comparison of models. 

We analyze the multiplicity distribution data first 
with single NBD function and then extend our s tudy 
with weighted superposition of two NBD functions [ij] , 
where single NBD resulted in poor agreement with 
data. Analyzing the data in terms of weighted superpo- 
sition of two distribution functions becomes pertinent 
because of the reported change in slope or appear- 
ance of sub-structure by the CMS experiment in multi- 
plicity distribution for n > 20 in its largest 77-interval, 
rjc < 2.4. According to the two-component model of 
Ref.- [llj], the shoulder- like structure in the multiplic- 
ity distribution of hadronic collisions at \/s= 0.9 and 
1.8 TeV could be explained by weighted superposition 
of two NBDs, representing two classes of events, "semi- 
hard - events with minijets or jets" and "soft - events 
without minijets or jets" . 

Analysis in terms of weighted superposition of two 
NBDs could be carried out following the formalism 
pH \2T\ of Giovannini and Ugoccioni for the LHC data, 
as has been suggested [131 by the authors. Beside 
proposing [Till the model to describe multiplicity distri- 
bution with sub-structure, Giovannini and Ugoccioni, 
predicted possible scenarios for structure of the "soft" 
and the "semihard" components of particle production 
in the full phase space in hadronic collisions in the 
TeV energy-range, in terms of the fc-parameters of the 
two-NBD and the clan parameters [28|. The study of 
multiplicity distribution and the clan structure analysis 
psl j in terms of the two-component model was then ex- 
tended [131 to the limited pseudorapidity intervals for 
the situation when analysis in limited phase space is 
carried out after classification of events into "soft" and 
"semihard", to ensure that the weight factors for the 
components in the limited phase space remain same as 
those in the full phase. The resulting function of the 
weighted superposition of two NBDs is given by: 

Pn[{n)soft{\fs, Vc),ksoftiVs,Vc)] + [1 ~ asoftiVs)] 

where a^o/t is the fraction of "soft" events and is a 
function of ^/s only. The other parameters, functions 
of both, the ^/s and the rjc, have usual meanings as 
described for Eq. - (1) with suffixes in parameters in- 
dicating respective components. 

In the discussed model, the parameter ksoft is con- 
stant with energy of hadronic collisions, indicating va- 
lidity of KNO-scaling [3l| for the "soft" component. 



In terms of ksernihard, toT the "semihard" component, 
which is expected to be dominant in the TeV energy 
range of LHC, the model proposes three scenarios: 

1) the "semihard" component also follows the KNO- 
scaling, where ksemihard remains constant with energy 

2) ksernihard dccreases linearly with increasing energy 
indicating violation of KNO-scaling by the "semihard" 
component and 3) a QCD-inspired scenario, where 
the KNO-violation is not as strong as in scenario 2, 
ksernihard starts decreasing with energy but asymptoti- 
cally tends to a constant value. The energy dependence 
of {n)soft and (n) sem.ihn.rd in the TeV energy range has 
been extrapolated [ll|] empirically from the data in the 
GeV energy domain. 

By fitting the LHC data with two-NBD function, 
constrained with the predictions for {n)soft and ksoft 
as in Ref. [l^] , the parameters related to the "semihard" 
component could be obtained to match with the re- 
spective predictions. But, as the data, available with 
us for the analysis, contains both the "soft" and the 
"semihard" components and there is no published LHC- 
data, as yet, to analyze the "soft" and "semihard" sub- 
samples separately, the parametrization of Giovannini 
and Ugoccioni cannot be used. 

In the context of the superposition of the "soft" and 
the "semihard" components, some significant observa- 
tions [H, H^, 1131 in collider experiments at SPS, Teva- 
tron and RHIC are worth mentioning. The analysis [l^ 
of data of pp collisions at ^/s= 630 and 1800 GeV by 
CDF experiment at Tevatron, Fermilab in two isolated 
sub-samples of "soft" and "hard" events revealed invari- 
ance of properties of "soft" sub-sample as a function of 
^/s. The energy invariance of dynamical mechanism of 
inelastic multiparticle production in "soft" pp collisions 
has been observed [30] to be valid at y/s~ 200 GeV by 
the STAR experiment at RHIC, BNL also. A compara- 
tive study 12] of charged particle multiplicities arising 
from non-single diffractive inelastic hadronic collisions 
at 30 GeV to 1800 GeV, including the data at 

cohider energies at ^5= 200 GeV to 1800 GeV in UA5 
(SPS) and E735 (Tevatron) experiments, revealed that 
the multiplicity distribution data of collider energies de- 
viate from much discussed KNO-scaling [sij , which sat- 
isfactorily explains multiplicity distributions up to the 
ISR energies. The deviation, in the form of a shoulder- 
like structure apparently appeared in the collider data 
due to superposition of distribution of particles from 
some other process, different from the KNO producing 
process, on the top of the KNO distribution. 

The "hard" events in the referred experimental anal- 
ysis and the "semihard" events as termed in the dis- 
cussed model are essentially similar class of events, in- 
volving hard parton-parton scatterings (due to high 
momentum transfer) resulting in QCD jets of high 
transverse momentum above a certain threshold. 

In terms of two-NBD, we aim to study the energy de- 
pendence of the two components of particle production. 
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Our interest lies particularly with the "soft" component 
of particle production in pp collisions at LHC, in view 
of the energy invariance of the "soft" component ob- 
served in collider energies prior to LHC. We extend the 
study further to the clan structure analysis [13, H^l for 
the two components. 



IV. ANALYSIS AND DISCUSSIONS 



A. Behavior of the NBD parameters 



^/i (TeV) 




{n)NBD 


(n) 


x'/d.o.f 


0.9 


0.5 


3.66 ± 0.04 




6.12/22 




1.0 


7.49 ± 0.08 


7.26«;li 


53.35/38 




1.5 


11.32 ±0.10 


10.95l»;f, 


49.55/50 




2.0 


15.26 ±0.13 


14.831";?^ 


36.69/60 




2.4 


18.36 ±0.14 


17.86lE;;fo 


46.29/66 


2.36 


0.5 


4.70 ± 0.08 




6.38/21 




1.0 


9.42 ±0.11 


9.26«:1? 


55.30/38 




1.5 


14.35 ±0.16 


14.011^2? 


24.79/48 




2.0 


19.35 ±0.21 


18.93l»i? 


29.11/58 




2.4 


23.35 ±0.25 


22.63l«;i 


29.76/68 
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0.5 


6.16 ±0.05 


5.98t;j:i^3 


83.36/39 




1.0 


12.49 ± 0.08 


12.181";- 


152.65/68 




1.5 


18.89 ±0.10 


18.53lE;;f3 


226.57/93 




2.0 


25.47 ±0.14 


25.101°;?;, 


208.56/113 




2.4 


30.90 ±0.16 


30.32l°;^t 


129.37/125 



TABLE I: Table showing {n)NBD parameter of NBD, (n) 
calculated from published multiplicity distribution and the 
/d.o.f. for different pseudorapidity intervals for ^/s = 0.9, 
2.36 and 7 TeV. 



We fit the multiplicity distribution data 0, HI] for 
pseudorapidity intervals rjc — 0.5, 1.0, 1.5, 2.0 and 2.4 
with a single NBD function and tabulate mean mul- 
tiplicity (n) calculated from the distribution data, the 
values of best fitted parameter, {njNBD, and the cor- 
responding values of x^/d-o.f in Table-1. The increase 
in average multiplicity, (n) with increasing energy and 
with the width of the symmetric pseudorapidity inter- 
vals around ijcm = in the mid-rapidity region is a fact 
well established by experiments. As the parameter (n) 
of NBD function, given by Eq.- (1), gives the average 
multiplicity, the closeness of its best fitted values with 
respective measured mean multiplicity and its behav- 
ioral pattern with respect to ^/s and jTyj-intervals, as 
shown in Table-1, is expected. The behavior of (n)- 
parameter of NBD with respect to ^/s and rjc for the 
analyzed data is more clearly depicted in Fig.l. 

The Fig. 2 shows the dependence of the NBD- 
parameter, k on the size of the pseudorapidity interval, 
rjc and on the energy, y/s of collisions at LHC. In the 
considered energy domain, for a given rapidity inter- 
val, k decreases with increasing center-of-mass energy 



of collision. This behavior is consistent with the ob- 
served energy dependence of k in the pre-LHC energy 
range. 
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FIG. 1: Parameter (n) from NBD for <0.5 to 2.4 for 
= 0.9, 2.36 and 7 TeV. The error-bars associated with the 
data-points are not visible as the corresponding magnitudes 
are smaller than the dimension of symbol-size in the plots. 
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FIG. 2: Parameter k of NBD for \ri\ <0.5 to 2.4 for 
0.9, 2.36 and 7 TeV. 



On the other hand, for a given ^/s, though the general 
trend in the behavior of k shows increase in k with in- 
crease in the size of the symmetric pseudorapidity win- 
dow around the center-of-mass pseudorapidity rjcm — 0, 
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deviations appear at r]c = 1.0 for the ^/s = 0.9 and 2.36 

TcV data as have been shown in the Fig. 2. The reasons 
for such deviations is not understood. Further, the rate 
of increase in k with increased size of the pseudorapid- 
ity interval decreases with increasing energy. 




n 

FIG. 3: Primary charged hadron multiplicity distributions 
for <0.5 to 2.4 for = 0.9 TeV. The solid lines drawn 
along the data-points correspond to respective fits of NBD. 
The error-bars include both the statistical and the system- 
atic uncertainties 
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n 

FIG. 4: Primary charged hadron multiplicity distributions 
for \'n\ <0.5 to 2.4 for ^ = 2.36 TcV. The sohd hncs drawn 
along the data-points correspond to respective fits of NBD. 
The error-bars include both the statistical and the system- 
atic uncertainties 



For \/s = 7 TeV, however, the trend is followed up to 
\ri\ <1.5 beyond which the trend gets reversed. The 
reason for this deviation may be due to the appearance 
of sub-structure in the multiplicity distributions and 
so the inadequacy of a single NBD function to fit the 
distributions as has already been reflected in terms of 
the x^/d-o.f values, as listed in Table - 1, for the fits 
to the distributions at ^/s = 7 TeV. 




n 



FIG. 5: Primary charged hadron multiplicity distributions 
for <0.5 to 2.4 for = 7 TeV. The solid lines drawn 
along the data-points correspond to respective fits of single 
NBD while the dashed lines correspond to respective fits of 
Two-NBD. The error-bars include both the statistical and 
the systematic uncertainties 



We, therefore, fit the ^/s = 7 TeV data with weighted 
superposition of two NBD functions (given by Eq. - (3)) 
also. We plot multiplicity distribution data for ^/s — 
0.9 TeV (Fig.3), for ^ = 2.36 TeV (Fig.4) along with 
best fitted NBD function and for = 7 TeV (Fig.5) 
along with best fitted NBD function and superposition 
of two NBD functions. 

As could be seen in Fig.- 5, the multiplicity distri- 
bution data of charged hadrons for ^/s= 7 TeV for all 
the pseudorapidity intervals fit better to the weighted 
superposition of two-NBDs than to a single NBD. The 
improvement is more clear in terms of ^ j d.o.j values 
which, along with values of best fitted free parameters 
of the two-NBD function given by Eq. - (3), are tabu- 
lated in Table - II. From Table-II, we also observe sys- 
tematic trend in 77c-dependence of the best fitted values 
of parameters of two-NBD, for ^/s = 7 TcV data. We 
find {n}semihard ~ Hn) soft for a given 77-interval. kgoft 
increases and ksemihard decreases with increase in r]c- It 
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{n)s 




{n)sh 




1.30 ±0.29 
1.48 ±0.28 
1.73 ±0.28 
1.98 ±0.27 
2.38 ±0.34 


4.53 ± 1.11 
8.90 ± 1.67 
12.16 ± 1.89 
14.86 ± 1.89 
15.06 ± 1.48 


5.61 ±4.17 

5.22 ± 1.17 
4.74 ± 0.87 

4.23 ±0.72 
3.25 ±0.49 


14.08 ±2.32 
26.37 ±3.26 
36.35 ± 3.89 
44.87 ±4.12 
46.86 ± 3.45 


4.14/35 
6.97/64 
12.92/89 
14.93/109 
11.91/121 



Interestingly, however, for the largest available rj- 
interval, where shoulder-like structures appear Q in 
multiplicity distributions, the unconstrained free pa- 
rameters of the two-NBD show systematic behavior 
with respect to ^/s, as can be seen in Table - III. 
The lack of systematic trend in two-NBD parameters 



TABLE II: Table showing best fitted parameters of function 
given by Eq. - (3) and corresponding x^/d-O.f. for multi- 
plicity distributions in pseudorapidity intervals rjc <0.5, 1.0, 
1.5, 2.0 and 2.4 (tabulated in order) for ^ = 7000 GeV. 
Suffixes s and sh to the title of columns represent soft and 
semihard, respectively. Soft events fraction decreases from 
84% {ric <0.5) to 49% (r?c <2.4). 



may be noted that the values of x^/d-o.f become lower 
than unity which may be due to the fitting of the func- 
tion considering data-points with large errors (include 
both statistical and systematic errors). 



B. Energy invariant multiplicity distribution of 
soft component 

As already discussed in Section - III, the study of 
multiplicity distributions in terms of two-NBD could be 
important in the context of the experimental observa- 
tions j29l . Isoj of energy invariance of multiplicity distri- 
butions of the soft component of events. We have also 
discussed why the predictions [l^ by Giovannini and 
Ugoccioni on two-NBD parameters for the TeV energy- 
domain in limited phase space could not be used for 
our analysis. Alternately, the two-NBD with parame- 
ters related to the "soft" component constrained with 
low-energy data could have been used. But, as the two- 
NBD parameters are functions of both the ^/s and the 
rjc, the study of dependence on ^/s could be possible 
only in the same rjc and vice versa. 

So, to study the ^/s and rjc dependencies of the two 
components, we prefer to study the distribution data 
of CMS at = 0.9, 2.36 and 7 TeV in 7]^ < 0.5, 1.0, 
1.5, 2.0 and 2.4 in terms of the two-NBD function given 
by Eq. - (3) with unconstrained parameters (keeping 
ksoft, {n)soft, ksermhard and (u) semihard all free to pro- 
duce best fit). All these fittings resulted in lowering 
of /d-o.f (to different extents) as compared to fits 
with single NBD, although the improvement in the fit 
by two-NBD is not that significant in small rjc at ^/s 
= 0.9 and 2.36 TeV. It is important to note that fit- 
ting the distributions as measured by ALICE ^ up to 
77c < 1.3 at ^/s = 0.9 and 2.36 TeV with sum of two 
NBDs also did not significantly improve the description 
of data. Also, the values of best fitted free parameters 
obtained from our two-NBD fits at small pseudorapid- 
ity intervals (rjc), do not follow any systematic trend 
with respect to ^/s or rjc, in general, for the data at ^/s 
= 0.9 and 2.36 TeV . 



y^{TeV) 




{n}s 


^sh 


{n)sh 


0.9 
2.36 
7.0 


2.44 ± 0.32 
2.57 ±0.52 
2.38 ±0.34 


14.78 ± 1.99 
15.74 ± 2.98 
15.06 ± 1.48 


8.13 ±2.34 
6.27 ±2.21 
3.25 ± 0.49 


35.11 ±3.90 
41.92 ±6.21 
46.86 ± 3.45 



TABLE III: Table showing best fitted free parameters of 
function given by Eq. - (3) for multiplicity distributions in 
pseudorapidity interval, r/c <2.4 for ^/s — 0.9, 2.36 and 7 
TeV. Suffixes s and sh to the title of columns represent soft 
and semihard, respectively. Soft events fraction is 82% for 
= 0.9 TeV, 69% for -/s = 2.36 TeV and 49% ^ = 7 
TeV. 



at small r]c in contrast to that at large r]c, where the 
energy invariance of parameters related to soft compo- 
nent is observed in the LHC data, needs further discus- 
sion, particularly, as the energy-invariance of the "soft" 
component of particle productions has been observed in 
77c <l-0 by both the STAR and the CDF experiments 
and the model involving superposition of two NBDs has 
been shown to be valid in small 77c-intervals. The dif- 
ferent behavior of LHC data at small 77c-interval may 
be due to analysis of different class of event samples as 
compared to STAR and CDF. We recollect, the valid- 
ity of the two component model of Giovannini and 
Ugoccioni for the full phase space was extended [27,] to 
limited pseudorapidity intervals only after classification 
of events is carried out in the full phase space, to ensure 
multiplicity distributions in limited intervals holding 
same weighting factor as that in the full phase space. 
Also, both the STAR and the CDF experiments ana- 
lyzed isolated "soft" and "hard" sub-samples of events 
separately. But, this work deals with available LHC 
data sample including both the "soft" and the "semi- 
hard" events for a given 77-interval and attempts to ex- 
tract contributions of the components by fitting mul- 
tiplicity distribution data by weighted superposition of 
two NBDs. A small pseudorapidity interval may not 
include all the particles of an event and possibility of 
exclusion of part of an event in the small interval is 
more for high multiplicity events, which are likely to 
be abundant at new LHC-energies. So, the multiplic- 
ity distribution of a sample of all events, including the 
"soft" and the "semihard" components, in small inter- 
val of pseudorapidity may not reflect the same weight- 
ing factor for the two-NBD fit as that in the full phase 
space. The effect of partial exclusion of an event min- 
imizes with increasing size of phase space. In a large 
pseudorapidity interval, therefore, the weighting factor 
of the two-NBD may remain same as that in the full 
phase space. 
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Nevertheless, application of the formalism of the 
weighted superposition of two NBDs in the analysis of 
the published multiplicity data of CMS reveals signifi- 
cant property of energy invariance of "soft" component 
of particle productions at LHC in the largest available 
pseudorapidity interval, rjc <2.4, where shoulder-like 
structure appear in the multiplicity distributions. At 
this point, it is important to compare the goodness of 
the fits to the multiplicity distribution data by a single 
NBD and by the superposition of two-NBDs. To com- 
pare the goodness of fits at the considered LHC energies 
in the pseudorapidity interval, 77c < 2.4, we carry out 
the residual analysis andplot the residuals in Fig. - 6. 
The residual is defined as the difference between a 
data-point and corresponding fit-value. It is clear from 
the plots in Fig. - 6, that the weighted superposition of 
two NBD functions fits better than a single NBD with 
the available multiplicity distribution data of LHC at 
= 0.9, 2.36 and 7 TeV at 77^ < 2.4. 




FIG. 6: Plots of residual analysis to test the goodness of 
fits for = 0.9 TeV, = 2.36 TeV and = 7 TeV. 
The dotted lines correspond to fits to a single NBD func- 
tion and the continuous lines correspond to fits to weighted 
superposition of two NBDs. 



Clan structure analysis in the two-component 
model 



Better agreement of two-NBD with data at y/s— 7 
TeV in all available rjc and at y/s = 0.9 and 2.36 TeV 
in < 2.4, motivates us in contrasting the LHC data 
at TeV energies with the clan structure analysis within 
the framework of the two-component model [ll], [23| • 

The discussed two-component model fits in the 
framework of clan structure analysis [2^. In fact, the 
authors of Ref. |27| have pointed out that to distin- 
guish different scenarios, one should look at the 1/k pa- 
rameters and related clan structure analysis. The clan 
model [m is based on cascading process, where par- 
ticles are emitted from a previously produced particle 
while the producing particle can change its momentum 
and quantum numbers during the process, as it hap- 
pens in case of well known fragmentation and decay 
processes. The group of particles including one origi- 
nally produced from the collision, directly or indirectly, 
and particles produced from that in following steps of 
cascading form a cluster with the originally produced 
first particle termed as the ancestor of the cluster. Such 
a cluster or a group of particles with common ancestry 
is termed as clan |28| . As per definition, a clan contains 
at least one particle. Clans can be assumed to be pro- 
duced independently. The characteristic parameters of 
the clan model are the average number of clans, N and 
the average number of charged particles per clan, fic, 
which are related to the NBD parameters as follows: 



and 



N ^kx ln(l + ^) 
k 



N 



(4) 



(5) 



For the two-NBD two-component model, describing 
"soft" and "semihard" components of events, the be- 
havior of the clan parameters need to be studied sepa- 
rately for each component. 

For the study of -yi-dependence of clan parameters, 
we consider the multiplicity distributions and corre- 
sponding fit-parameters of two-NBD in 77c < 2.4 at y/s 
= 0.9, 2.36 and 7 TeV and plot the clan parameters, 
along with the respective NBD-parameter 1/k, for both 
the "soft" and the "semihard" components in Fig. - 7. 
For the study of ?7c- dependence of the same parameters, 
we consider multiplicity data in r/c < 0.5, 1.0, 1.5, 2.0 
and 2.4 at y/s= 7 TeV and plot, in Fig. - 8, the param- 
eters similar to those in Fig. - 7. The NBD-parameter 
1/fc is plotted along with the clan parameters for con- 
venience of comparing the behavior of all the three pa- 
rameters, 1/fc, 77c and the N with predictions of the 
two-component model [III [l^l . 

As can be seen from the plots in Fig. - 7, for the 
"soft" component, all the three parameters, character- 
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izing different scenarios as particle productions accord- 
ing to the discussed model in Ref.-jl^lj the NBD pa- 
rameter, the average number of charged particles 
per clan, fic and the average number of clans, N are 
constant with energy. For the "semihard" component, 
1/k and fic increase while N decreases with increase 
in energy. The rate of change in the parameters is 
rapid in the range from 2.36 TeV to 7 TeV as com- 
pared to the that in the range from 0.9 TeV to 2.36 
TeV. Comparing the nature of the energy-dependence 
of these parameters (though not for same r]c) as shown 
in Ref. [23| , one sees that the behavior of energy depen- 
dence of clan parameters and the fcsemihard-parameter 
of two-NBD in 77c< 2.4 in the new LHC-energies match 
with the scenario - 2 (as discussed in Section - III) of 
the two-component model of particle production, repre- 
sented by weighted superposition of two NBD, indicat- 
ing violation of KNO scaling by the "semihard" com- 
ponent. But, contrary to the prediction by the model, 
^/ksoft is always larger than 1/ksemihard in the consid- 
ered energy-range. 
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FIG. 7: Energy dependence of the NBD parameter 1/k 
(panel a), the average number of charged particles per clan, 
fic (panel b) and the average number of clans, N (panel c) 
obtained for "soft" and "semihard" components from mul- 
tiplicity distribution data in |ry| <2.4 in the LHC energy- 
domain {^/s = 0.9, 2.36 and 7 TeV). The lines in the plots 
are drawn joining the points to guide the eye. 



The behavior of ?7c-dependence, as depicted in Fig. - 
8, of the characteristic parameters for the two com- 



ponents, however, do not corroborate the finding from 
the study of the energy dependence. According to the 
model, the f /kgoft decreases and I / ksemihard decreases 
rapidly with increasing rjc in the scenario - 2 of the 
model. The plot of rjc - dependence of l/ksemihard for 
the data at ^/s shows rising tendency. 
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FIG. 8: The 77c-dependence of the NBD parameter 1/k 
(panel a) , the average number of charged particles per clan, 
fic (panel b) and the average number of clans, iV (panel c) 
obtained for "soft" and "semihard" components from mul- 
tiplicity distribution data at y/s — 7 TeV. The lines in the 
plots are drawn joining the points to guide the eye. 



V. SUMMARY 

Our study on multiplicity distributions of charged 
hadrons from pp collisions in limited pseudorapidity in- 
tervals at LHC energies in terms of Negative Binomial 
Distribution function reveals the followings: 

f ) At pseudorapidity intervals of small sizes, particu- 
larly at energies ^/s = 0.9 and 2.36 TeV, a single NBD 
function fits the distribution data reasonably well, while 
parameters of two-NBD show no systematic trend. For 
the distribution data at ^/s = 7 TeV, however, a sin- 
gle NBD function appears inadequate, while weighted 
superposition of two NBDs fit the data satisfactorily. 

2) The energy {^/s) invariance of the parameters re- 
lated to the "soft" component and so of the respective 



9 



multiplicity distribution at rj^ < 2.4, where the mea- 
sured distributions show shoulder-like structure. This 
observation could be indicative to the invariance of the 
dynamical mechanism of "soft" multiparticle produc- 
tion, as has already been seen by the STAR experiment 
in pp collision data at -^/s — 200 GeV and by the CDF 
experiment in pp collisions at y/s = 630 and 1800 GeV. 

3) The multiplicity distributions for all the available 
LHC energies in 77c < 2.4 agree better with weighted 
superposition of two NBDs than a single NBD function. 

4) Behavior of clan parameters show energy- 
dependence in accordance with one of the predicted 
scenarios by the two-source model. But the study of 
77c-dependence of data at y/s = 7 TeV, which otherwise 
fit with two-NBD, does not substantiate the finding in 
the energy-dependence study. 

This study of multiplicity distributions of the pp col- 
lisions at LHC, in terms of NBD, highlights significant 
features of mutiparticle productions in hadronic colli- 
sions at LHC energies. In spite of limitations in the 
available data in small pseudorapidity intervals vis-a- 
vis the adopted two-NBD formalism, the two-NBD de- 
scribes the multiplicity distributions data better in all 
pseudorapidity intervals at \/s~ 7 TeV and in large 
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